Augmented Recursion For One-loop Amplitudes 



David C. Dunbar, James H. Ettle and Warren B. Perkins 
Department of Physics, 
Swansea University, 
Swansea, SA2 8PP, UK 

We present a semi-recursive method for calculating the rational parts of one-loop amplitudes when recursion 
produces double poles. We illustrate this with the graviton scattering amplitude M 1_loop (l _ , 2 + , 3 + , 4 + , 5 + ). 



1. Introduction 

On-shell recursive techniques, using the ratio- 
nality of tree amplitudes and their complex fac- 
torisation properties, have proven very success- 
ful in the computation of scattering amplitudes 
in gauge and gravity theories [TJ[2]. Specifically, 
in a theory with massless states, if we use a 
spinor helicity representation for the polarisation 
vectors it is possible to write the amplitude en- 
tirely in terms of spinorial variables A(X t a ,X\) 
where the massless momentum of the i th parti- 
cle is A^,A^ = (<Jfj,) a ak^. The analytic structure 
of the amplitude can be probed by choosing a 
pair a, b of external momenta and shifting these 
according to 

A a — ^X a -zX\ \ b — >X b + zX a (1) 

where we suppress the spinor indices. If the 
shifted amplitude A(z) (a) is a rational function, 
(b) has finite order poles at points zt, and (c) van- 
ishes as z — > oo, then applying Cauchy's theo- 
rem to A(z)/z with a contour at infinity yields 

poles Zi 

At tree level the factorisation of amplitudes 
is simple: amplitudes must factorise on multi- 
particle and collinear poles into the product of 
two tree amplitudes defined at z = Zi. Thus we 
can express the n-point tree amplitude in terms 
of lower point amplitudes, 

4r(o) = E<Ti CT ^)^2^-r+i(^). ( 3 ) 



where the summation over i is only over factori- 
sations where the a and b legs are on opposite 
sides of the pole. This technique is very effec- 
tive in computing tree amplitudes and has been 
extended to a variety of other applications includ- 
ing gravity [2]. 

Beyond tree level there are three potential bar- 
riers to using recursion. Firstly, the amplitudes 
generally contain non-rational functions such as 
logarithms and dilogarithms; secondly, the am- 
plitudes may contain higher-order poles for com- 
plex momenta; and finally, the amplitudes may 
not vanish asymptotically with z. Nonetheless a 
variety of techniques based upon recursion and 
unitarity have been developed. A one-loop am- 
plitude for massless particles may be expressed 
as 

A i-ioo P= j2 cai + R, (4) 

n=2,3,4;i 

where the scalar integral functions l l n are the var- 
ious scalar box, triangle and bubble functions. 
The amplitude can thus be determined by com- 
puting the rational coefficients, Ci, and the purely 
rational term R. The c 2 ; can be computed by the 
four-dimensional unitarity technique [3H5] or in- 
deed recursively [6|. Many techniques have been 
developed for evaluating R: D-dimcnsional uni- 
tarity, recursion and specialised Feynman dia- 
gram techniques [7H2"U]. 

In general, the rational term R does not sim- 
ply satisfy the requisites for recursion. If the am- 
plitude has only simple poles but does not van- 
ish as z — > oo then it can be possible to for- 
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mulate auxiliary recursion relations [21]. How- 
ever there are rational amplitudes for which one 
cannot find a shift which only generates sim- 
ple poles such as the single-minus amplitudes 
A x - loo P(l-, 2+, • • • ,n+). These amplitudes van- 
ish at tree level and consequently are purely ra- 
tional at one-loop. A shift on these amplitudes 
yields double and single poles. The double pole 
is not in itself a a barrier to using recursion, how- 
ever to obtain the full residue one needs to know 
the coincident single pole, or the 'pole under the 
double pole', which is not determined by factori- 
sation into on-shell amplitudes. In [22] for Yang- 
Mills this was postulated to be 



S(a 1 ,K+,a 2 )S(b 1 ,K-,b 2 ) 



(K 2 ) 2 



K 2 



(5) 



S(a, s + ,b) = 



where the 'soft' factors arc 
(ab) I ((a s) (sb)), S(a,s~,b) — [a b] / ([a s] [a b}) , 
and ax, a 2 (£>i, b 2 ) are colour-adjacent to K on 
the left (right) side of the pole. With this ansatz, 
recursion correctly reproduces the known single- 
minus one-loop amplitudes. In [23 it was shown 
that the consistency requirements for recursion in 
QCD arc sufficient to determine these soft factors. 

The above postulate, or variations thereof, does 
not work for gravity amplitudes [21] ■ Here, 
we apply a semi-recursive technique for grav- 
ity scattering amplitudes that obtains the 'pole 
under the pole' using an axial gauge formalism 
to calculate the previously-unknown amplitude 
M 1 - 1o °p(1-,2+,3 + ,4+,5 + ). We assume that the 
shifted amplitudes vanish as z — > oo. The 
derived amplitude has the correct symmetries 
and soft limits, providing strong evidence for the 
validity of this assumption. Further, we have 
checked the result by a completely independent 
'string-based rules' [25ll2fj] computation. 

2. Recursion 

The factorisation of one-loop massless ampli- 
tudes is described in [27] . 



i l-loop 
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l-loop 



r+1 
l-loop 



K 2 
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n— r+1 



4 tree - A 
A r+1 K 2 ^-n-r+l 



A tree 4 tree p 

K 2 n_, '+ 1 



(6) 



where the one-loop 'factorisation function' F n is 
hclicity-independent. Naively this only contains 
single poles, however for complex momenta there 
are double poles. These can be interpreted as due 
to the three-point all-plus (or all-minus) one-loop 
amplitude also containing a pole 



* l-loop 
A 3 



M) = ±- 2 V^° op {K+,a\b+) 



(J) 



where, for pure Yang-Mills, 



V. 



l-loop 



(K+,a+,b + ) = - 



487T' 



[Ka] [ab][bK]. 



(8) 



Explicitly, consider the amplitude 



i l-loop 



(l-,2+ 3+4+ 5+) 



(34Y 



[2 5] 



[12] [51] 



(14) 3 [4 5] (3 5) (13) 3 [3 2] (4 2) 



(12) (2 3) (4 5}' <15><5 4><3 2>' 



(9) 



If we carry out a complex shift on As,Ai as in 
eq. ([!]) then (4 5) — > (4 5) + z (4 1) which van- 
ishes at z — —(4 5)/ (41) and the amplitude has 
a double pole at this point. 

Computing this amplitude using V^~ loop cor- 
rectly generates the double pole in the ampli- 
tude 22 ): 24], however it needs augmentation to 
give an expression with the correct single pole. 
By trial and error, adding the second term in ([S]) 
gives the correct single pole and completes the 
computation of the amplitude. 

For gravity the vertex 



V 1 - loop (K + ,a + ,b+) 



iK 3 {[Ka] [ab] [6 if]) 2 



1440tt 2 



(10) 



can be used to generate a double pole term but 
attempts [23] to implement a universal correction 
for the single pole analogous to that of ([5]) have 
failed. The resolution is to replace the factori- 
sation term of © with a tree insertion diagram: 
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which we compute using axial gauge diagrammat- 
ics. The circle in the diagram represents the sums 
of all possible tree diagrams with two internal legs 
and the given external legs, which we denote r. 
Note that we evaluate these diagrams for real mo- 
menta and only carry out analytic shifts on the 
final expressions. 

3. Axial gauge diagrammatics 

Following 29 a we use a set of Feynman rules 
for Yang-Mills amplitudes based on scalar prop- 
agators connecting three and four point vertices. 
The starting point is the expansion of the axial 
gauge propagator in terms of polarisation vectors, 



When adopting a recursive approach which 
involves shifting a negative-helicity leg a and 
a positive-helicity leg b, the recursion-optimised 
choice for the reference momentum q is 



A„ — A?,. 



(15) 



With this choice of q the leg a (6) can only enter 
a diagram on a V 3 MHV (Kj MHV ) vertex, and there 
are no four-point vertices in the single-minus am- 
plitudes at tree or one-loop level. 

Singularities arise in the loop integration from 
the region of loop momentum where the denomi- 
nators of three adjacent propagators vanish simul- 
taneously. This requires the two null legs to which 
the propagators connect to become collinear. In 
the integration region of interest all the legs of 
r are close to null and t approaches the corre- 
sponding on-shell tree amplitude. The internal 
legs are also close to collinear. Helicity config- 
urations for which r is singular in this collinear 
limit, shown in Fig. 1, contribute to the double 
(and single) pole, conversely those that give a 
vanishing r in the collinear limit give no residue. 



(11) 



where 



V2~(k»q)' M V2[k>>q] 



2k ■ q 



with 

k b := k 



2k ■ q 



(12) 



(13) 



where q is a null reference momentum which may 
be complex. The resulting three-point vertices 
are, 

n 2 



1 

i'a/2 



V, 



MHV/i- o- 



(l-,2-,3+) 



(12) [3qY 
[U] [2 5] 
[21] (3g) 2 
(1?) (2?) 



(14) 



_^y 3 MHV (1+)2+)3 - ) 
i\l 2 

along with a V3(l + ,2~,3°) vertex which can be 
absorded into effective four-point vertices. 





(a) (b) 
Figure 1: Contributing helicity structures 

The diagram of Fig. 1(b) evaluates to 



dV 



[b|/|a)[c|f|a) (Co) 2 r(C + ,d+,--- ,o~,B~) 
(ba)(ca) ( Ba ) 2 l 2 (l + k b ) 2 (l - k c f 

(16) 



where B = I + b, C — c — I, and the momenta in 
the spinor products are q- nullified as in (fl3|) . We 
construct a basis for the loop momentum using b 
and c: 

I = ai(k b + fc c ) + a 2 (k b - k c )+ 

(a 3 + tai)— — AfcA c + (a 3 - ion)- — rA c A 6 
{ba) {caj 
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Under this parametrisation, 



dH f(l) 



i 2 (i + k b y{i-k c y Sbc 



da, F(ai)f(l(ai)) 
(18) 



where F(cti) has no dependence on Sf, c . The in- 
tegrand from Fig. 1(b) then becomes, 



[be] (CaY 
(be) (BaY 



T(C+,d+,...,a-,B-) x F'(ai). (19) 



In order to evaluate the contribution from (|19p 
we must evaluate the tree structures to order 
(bey. For diagrams within r involving l/sf, c , this 
means going beyond leading order. These corre- 
spond to triangles in the full diagram and the 
calculation is readily done exactly. The diagrams 
without this propagator need only be calculated 
to leading order. In this regard, not only is the 
recursive approach selecting a subset of diagrams 
for calculation, it is also allowing us to calculate 
these diagrams in a very convenient limit. 

For gravity the equivalent expression to (fT!?)) is 



[be] 3 (Co) 4 
(be) IB a 



^T s (C + ,d + ,...,a-,B-)xF( ai ). (20) 



4. The graviton scattering 

M 1 - loop (l-,2+,3 + ,4+,5 + ) 



amplitude 



There are three types of recursive contribution 
to this amplitude, which in turn are summed over 
the distinct permutations of c, d and e. Dia- 
grams R\ and R2 involve only single poles and 
are obtained from the corresponding four-point 
one-loop amplitudes for the circles marked L. 
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Doing recursion with the shift (JXJ) , we obtain 

_ 1 (ad) 2 (ae) 2 [bc][de} 4 
l(a ' ' C ' ' 6) 5760 (ab) 2 (bc)(ce) 2 (cd)i(de) 2 X 
((cd) 2 (ae) 2 + (ac)(cd)(de)(ae) + (ac) 2 (de) 2 ) , 

(21) 

3 (ae)[be]±_ x 

(22) 



R 2 (a,b, c,d,e) 



5760 (cd) 2 [ab} 2 [ae] 
([bc] 2 [de] 2 + [bc}[cd}[de}[be] + [cd] 2 [be] 2 ) 



Diagram R3 contains a double pole so we must 
evaluate r g of (|20|) . We use the five-point KLT 
relation [50] . 

M(a-B-C + d+e+) = 

SBcSdeA(a-B-C + d+e + )A(a-C + B-e + d + ) + 
s Bd s C eA(a-B-d + C + e + )A(a-d + B-e + C + ) 

(23) 

in a form that restricts the (b c) pole to the first 
term. We calculate this as Laurent series in (b c) , 
dropping terms that will not contribute to the 
residues. While the KLT relations are only valid 
for on-shell momenta, we assume the deviation 
of ([2"3")l from a direct off-shell calculation may be 
neglected in the region around B 2 = C 2 = 0. 

In our choice of axial gauge, A(a~ B~ C + d + e + ) 
receives contributions from five diagrams, only 
two of which contain a Vs(B^ , C + , x) vertex and 
thus contribute to r's collinear singularity. De- 



1 The general case is worthy of further study I31| . 



5 



noting these by D a and Db, we find, using (|T7)) 



D n + D h 



(a\bc\a) 



(Co) 2 Sbc[ab](da) (ea) 
'[b\ad\e] - [6|cfe|e] 



[ae] (de) 



/«(«*), (24) 



where f a (oti) is some function that depends only 
on the integral parameters, a*. We note that the 
second term is sub-leading in the (6 c) pole. 

The leading pole in A(a~C + B~e + d + ) is ob- 
tained similarly and we obtain the full leading 
(be) pole in ([23]) as 



(Ba) 4 (ab) 2 (ac) 2 [def[bc] . . 

-SbcSde „ \ / t . \ r . i .7 , rTJa( a i)- 



(Ca)< 



(bc)(de)[a\d- 



(25) 



Combining this with the factors arising from the 
left hand part of the full diagram (c/. ([20]) ) and 
integrating over the on the leading term in the 
Laurent series for R% is proportional to 



[bc] 4 (ab) 2 (ac) 2 [de] 3 
(bc) 2 (de)[a\d + e\a) 



(26) 



which clearly displays the double pole factor. 

We now express each sub-leading contribution 
to ([2U)l as V x Sjfj(ai). Firstly there is the sub- 
leading contribution of ((24)) together with the cor- 
responding contribution from A(a~ C + B~ e + d + ): 



Si = 



s bc [be] s bc [bd] 



(27) 



[b\ad\e] [b\ae\d]' 

The remaining diagrams for A(a~ B~C + d + e + ) 
(and its counterpart A(a~C + B~e + d + )), in which 
B and C enter on different vertices contribute 



s bc [e\a\c) 
s a b[e\d\c) 



(28) 



^ = (bc)(de) f [e[B[a)[eb] | \d\B\a) [db] \ 
s ab [de] \ (da)(cd) (ea)(ce) J' 

These diagrams are finite in the collinear limit, 
so we can drop terms proportional to B 2 and C 2 . 
Finally we need the second term in (|2"5|) , which is 
also finite in the collinear limit and can be evalu- 
ated using MHV tree amplitudes, yielding: 



^4 = 



(bc)(de)[d\B\a)[e\C\a) 
[bc][de}(ab) 2 (cd)(ce) 



(30) 



Up to O ((be) 1 ) (|2H|) is then expressed as 

[bc} 4 (ab) 2 (ac) 2 [de] 3 
(bc) 2 (de)[a\d + e\a) 



(31) 



This has purely polynomial dependence on the a^. 
The integration thus gives constant numerical fac- 
tors which may be obtained by direct evaluation 
or, more conveniently, by considering collinear 
limits. 

We now determine the amplitude recursively by 
applying the shift to the integrated (|3"Tj) and 
evaluating the residue at z = —(bc)/(ac). The 
coefficient of the double pole has a z dependence 
under this shift which generates a further contri- 
bution to the single pole since 



Res 



1 # 

Zi dz 



z=Zi z(z — Zi) 2 Z 2 

The full contribution from R3 is then 

1 (ab) 2 (ac) 4 [bc] 4 [de\ 



(32) 



i?3(a, b, c, d, e) 



where 



A(a, b, c, d, e) 



5760 (ad)(ae){bc) 2 {cd){ce){de) 
x (1 + A(a,b,c,d,e)), (33) 



I (ad) (be) l(ae)(bc) 



2(ab)(cd) 2(ab)(ce) 



_ \db][eb](bc)(de) _ [dc][ec](bc)(de)(ca) 2 
(dc) (ec) [be] [de] (dc) (ec) [be] [de] (ba) 2 
7 [dc] [eb] (be) (de) (ca) 7 [db] [ec] (be) (de) (ca) 



2 (dc) (ec) [be] [de] (ba) 2 (dc) (ec) [be] [de] (ba) ' 

(34) 

The full amplitude is the sum over contribu- 
tions arising from three orderings of external legs, 

M i-loo P ( 1 - ) 2 + t 3+, 4+, 5+) = R(l, 2, 3, 4, 5) + 

#(1,2,4,5, 3) + #(1,2,5, 3, 4), (35) 

(the full amplitude has a factor of zk 5 /167t 2 ), and 
each R is the sum of the recursive diagrams, 



R = R\ + #2 + R3 



(36) 
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This expresion has the correct collinear limits, is 
symmetric under interchange of pairs of positive- 
hclicity legs and agrees numerically with that cal- 
culated by string-based rules. We have also cal- 
culated M 1 - 1o °p(1-,2+,3+,4+ ! 5+,6+) 32 , and 
again checked that it has the correct symmetries 
and collinear limits. Mathematica code for the 
five- and six-point amplitudes may be found at 
http : //pyweb . swan. ac . uk/~dunbar/graviton. html 



5. Conclusions and remarks 

We have demonstrated how to augment recur- 
sion to determine the rational terms in ampli- 
tudes with double poles under a complex shift. 
Double poles are unavoidable in the case of the 
amplitudes A 1 - loo P(l-, 2+, 3+, . . . in both 

Yang-Mills and gravity. In the absence of a uni- 
versal soft factor analogous to ([3]) , to perform the 
augmented recursion the sub-leading poles must 
be determined on a case-by-case basis. While we 
have done this for both the five- and six-point 
single-minus gravity amplitudes, this procedure 
could be used to calculate any higher-point single- 
minus amplitude. 
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